#ifndef __TEST_DGBMV_H__
#define __TEST_DGBMV_H__

void testfunc(double matrix[],double matrix2[], double vector[], double loesung[]);

#endif
for(i = 0;i<1;i++){//time stepping
    if (i >0 ){
    dcopy_(&iter_z,H,&incx,H0,&incy);
    }

    for(j = 0;j<1;j++){
    ///solve function
    dcopy_(&iter_z,H,&incx,Hold,&incy);
    createFDmatrix(iter_z,dt,dz,material,phi,k0,H,H0,rhomatrix,AB,b,lda);
    dgbsv_(&iter_z,&kl,&ku,&nhrs,AB,&lda,ipiv,b,&iter_z,&info);
    if (info == 0){
        printf(" yes yes yes \n\n\n\n");
    }
    //cgmethod(A,b,H,iter_z,iter_z,lda);
    //solvefunc(iter_z,dt,dz,material,phi,k0,H,H0,solvecase,rhomatrix);
    dcopy_(&iter_z,b,&incx,H1,&incy);
    double alpha = -1;
    daxpy_(&iter_z,&alpha,Hold,&incx,b,&incy);
    value_norm = dnrm2_(&iter_z,b,&incx);
    printf(" picar %d \n ",j);
    if (value_norm < tol){
        dcopy_(&iter_z,H1,&incx,H,&incy);

        break;
    }else{
         dcopy_(&iter_z,H1,&incx,H,&incy);

    }
}
int s;
for(s = 0;s<iter_z;s++){
printf("Lösung%f ",H[s] );
}

}
